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ABSTRACT 

We present measurements of higher order clustering of galaxies from the latest 
release of the Canada-France-Hawaii- Telescope Legacy Survey (CFHTLS) Wide. We 
construct a volume-limited sample of galaxies that contains more than one million 
galaxies in the redshift range 0.2 < z < 1 distributed over the four independent 
fields of the CFHTLS. We use a counts in cells technique to measure the variance £2 
and the hierarchical moments S n — £ n /£2 (3 ^ n ^ 5) as a function of redshift 
and angular scale. The robustness of our measurements if thoroughly tested, and the 
field-to-field scatter is in very good agreement with analytical predictions. At small 
scales, corresponding to the highly non-linear regime, we find a suggestion that the 
hierarchical moments increase with redshift. At large scales, corresponding to the 
weakly non-linear regime, measurements are fully consistent with perturbation theory 
predictions for standard ACDM cosmology with a simple linear bias. 
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1 INTRODUCTION 

In our current paradigm, structures in the Universe origi- 
nate from tiny density fluctuations emerging from a primor- 
dial Gaussian field after an early inflationary period. Grav- 
itational instability in an homogeneous expanding Universe 
dominated by Cold Dark Matter indeed gives rise to the rich 
variety of structures seen today by hierarchical clustering. 
This picture has shown to be consistent with measurements 
of in million-galaxy surveys of the local Uni verse like the 
Sloan Digital Sky Survey (SDSSfl and 2dF (|Colless et all 



* Based on observations obtained with MegaPrime/MegaCam, a 
joint project of CFHT and CEA/IRFU, at the Canada-France- 
Hawaii Telescope (CFHT) which is operated by the National Re- 
search Council (NRC) of Canada, the Institut National des Sci- 
ence de l'Univers of the Centre National de la Recherche Scien- 
tifiquc (CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at Terapix available 
at the Canadian Astronomy Data Centre as part of the Canada- 
France-Hawaii Telescope Legacy Survey, a collaborative project 
of NRC and CNRS. 
f E-mail: wolk@iap.fr (MW) 
1 http://www.sdss.org/ 



l200ll ); however, our knowledge of the evolution of the dis- 
tribution of galaxies and its changing relationship with un- 
derlying dark matter density field remains incomplete. 

On small scales, the evolution of the galaxy distribution 
strongly depends on the physical processes of galaxy forma- 
tion, while on larger scales, both the initial power spectrum 
as well as to the global geometry of the Universe play an 
important role. Precise measurements of the galaxy distri- 
bution as a function of redshift and galaxy type therefore 
represents a fundamental tool in cosmology and the n-point 
correlation functions provide a particularly powerful method 
to characterise it. 

The two-point correlation function (along with is 
Fourier transform t he power spect r um) is the most widely 
used statistic (e.g. iPeeblesI ll98Ct iBaumgart fc Frvl Il99li ; 
iMartmezI 12009, for a general review on the subject) as it 
provides the most basic measure of galaxy clustering - the 
excess in the number of pairs of galaxies compared to a 
random distribution as function of angular scale. Never- 
theless, the two-point correlation function represents a full 
description only in the case of a Gaussian distribution for 
which all higher-order connected moments vanish by defi- 
nition. It is well kno wn that the actual galaxy distribution 
is non-Gaussian (e.g. lFrv fc Peebledll978hlSharp et al.lll984l : 
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ISzapudi et all 1 19921 ; iBouchet et ai1ll993l ; lGaztanagalll994h . 

Such a non-Gaussianity is already induced by nonlinear 
gravitational amplification of mass fluctuatio ns, even if they 



originated from an initial Gaussian field (e.g. Peebles! 
[pYvf Il984l ; ljuszkiewicz et al J 1 19931 ; IBernardeau et al 



1980 



2002 



and references therein). But it can also be present in initial 
conditions thereselves, f rom topological defects s uch as tex- 
tures or cosmic strings (Kai ser fc Stebbmslll984) or fluctu- 



ations coming from certain inflationary models ( Guth|[l98ll ; 
iBartolo et al.ll20o3 . and references therein). In addition, the 
process of galaxy formation results in a "biased" distribu- 
tion of luminous matter with respect to the underlying dark 
matter, which adds another level of complexity in the non 
Gaussian nature of t he process (e.g. iPolitzer fc Wiselll984l; 
Bardeen et all Il986l; fFry fc Gaztanagal Il993l ; iMo fc White] 
19961 ; IBernardeau et all |2002| . and references therein) . For 



all these reasons, high order moments are required to obtain 
a full statistical description of the galaxy distribution. 

In this paper we will use the method of counts in cells to 
investigate the evolution of the variance £2 and the hierarchi- 
cal moments S n = in/Q _1 (3 ^ n ^ 5) of the galaxy distri- 
bution from z ~ 1 to present day. Until now, the majority of 
papers which have looked at these quant ities have analysed 
large local-Universe reds hift surveys; s e elRoss et al.1 d2006al. 
l2007f) for th e SDSS and iBaugh et all (|2004T ); ICroton et al l 
l |2004l2007f ) for the 2dFGRS. Furthermore, only a few mea- 
surements have been made in the literature of the red- 
shift evolution of th e hierarchical moments. For example, 
ISzapudi et all (|200lh used a purely magnitude-limited cata- 
logue of ~ 710 000 galaxies with Iab < 24 in a contiguous 
4 x 4deg 2 region. By assuming a redshift distribution for 
their sources, they found that a simple model for the time 
evolution of the bias was consistent with their measurements 
of the evolution of £3 as a function of redshift. They con- 
cluded that these observations favoured models with Gaus- 
sian initial conditions and a small amount of biasing, which 
increases slowly with redshift. Using spectroscopic r edshifts 
from the VVDS-Deep survey, iMarinoni et al.l l|2008l ) recon- 
structed the three-dimensional fluctuation to z ~ 1.5 and 
measured the evolution of £ and S3 over the redshift redshift 
range 0.7 < 2 < 1. They found that the redshift evolution 
in this interval is consistent with predictions of first- and 
second-order perturbation theory. 

Measuring accurately the higher-order moments of the 
galaxy distribution requires both very large numbers of ob- 
jects and also a precise control of systematic errors and re- 
liable photometric calibration: no previous surveys have the 
required combination of depth and areal coverage to make 
reliable measurements from low redshift to z ~ 1 in large 
volume-limited sample of galaxies. 

In this work, we use data from the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS) which has a 
unique combination of depth, area and image quality, in ad- 
dition to an excellent and highly precise photometric cali- 
bration. 

The paper is organized as follows. In Sect. 2, we present 
our data set and the sample selection. In Sect. 3, we out- 
line the techniques we use to estimate the galaxy clustering 
starting from the two-point correlation function and going 
to higher orders. The results are described and interpreted 
in Sect. 4. Conclusions are presented in Sect. 5. 

Throughout the paper we use a flat ACDM cosmology 



(fi m = 0.23, n A = 0.77, H = lOOhkms- 1 Mpc -1 and cr 8 = 
0.8) with h = 0.71. All magnitudes are given in the AB 
system. 



2 OBSERVATIONS, REDUCTIONS AND 
CATALOGUE PREPARATION 

2.1 The Canada- France-Hawaii Telescope Legacy 
Survey 

In this work, we use the seventh and final version 
of the Canada-France-Hawaii- Telescope Legacy survey 
(CFHTLSfl The T0007 data release results in a major im- 
provement of the absolute and internal photometric cali- 
bration of the CFHTLS. Since T0007 is a reprocessing of 
previous releases, the actual area covered by this release is 
almost identical to previous versions (some additional ob- 
servations were added to fill gaps in the survey created by 
malfunctioning detectors). 

The improved photometric calibration has been 
achieved by applying recipes adopted by the Supernova 
Legacy Survey (SNLS) (SNLSQ for the SNLS/Deep fields to 
bot h the Deep and Wide fields of the CFHTLS (for details, 
see lRegnault et al1 l2009). This improved photometric accu- 
racy translates into improved photometric redshift precision 
and in a consistent photometric redshift calibration over the 
four disconnected patches of the CFHTLS Wide. 

Complete documentati on of t he CFHTLS-T0007 re- 
lease can be found at the CFHT0 site. In summary, the 
CFHTLS- Wide is a five-band survey of intermediate depth. 
It consists of 171 MegaCam deep pointings (of 1 deg 2 each) 
which, as a consequence of overlaps, consists of a total of 
~ 155 deg 2 in four independent contiguous patches, reach- 
ing a 80% completeness limit in AB of u* = 25.2, g — 25.5, 
r = 25.0, i — 24.8, z = 23.9 for point sources. The pho- 
tometric catalogs we use here were constructed by running 
SExtractoiQ in "dual-image" mode on each CFHTLS tile, 
using th e qri images to cr eate a x 2 - sc iuared detection im- 
age fsee lSzalav et al.lli999f ). The final merged catalogue is 
constructed by keeping objects with the best signal-to-noise 
from the multiple objects in the overlapping regions. Ob- 
jects are_selected using i-band MAG_AUT0 total magnitudes 
iKronl l|l980h and object colours are measured using 3 arcsec 
aperture magnitudes. 

After applying a magnitude cut at i < 22.5, the entire 
sample contains 2,710,739 galaxies. After masking, the fields 
have an effective area of 



S eS ,i = 58.33, 18.50, 37.04, 18.50 deg 2 



(1) 



for fields Wl, W2, W3, and W4 respectively, leading to a 
total area of 133 deg 2 . 



2 http : //www . cf ht . hawaii . edu/Science/CFHLS/ 

3 http://cfht.hawaii.edu/SNLS/ 



4 http : //www . cf ht . hawaii . edu/Sc ience/CFHLS/T0007/ 

5 http : //www. astromatic .net/sof tware/sextractor 
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2.2 Photometric redshift and absolute magnitude 
estimation 

The photometric redshifts we use here have been provided 
as part of the CFHTLS-T0007 release. They a re fully do<> 
umented on the corresponding pages at the TERAPIXj 
site and the photometric redshift release documentatior[j 
describes in details how the photometric redshifts were 
comp uted following |llbert et alj|2006l ; ICoupon et al]|2009t 
120121 ) . 

The photometric redshifts are computed using the 
"Le Phare'{f| photometric redshift code which uses a stan- 
dard template fitting procedure. These templates are red- 
shifted and integrated through the transmission curves. The 
opacity of the intra-galactic medium (IGM) is accounted for 
and internal extinction can be added as a free parameter 
to each galaxy. The photometric redshifts are derived by 
comparing the modelled fluxes and the observed fluxes with 
a x 2 merit function. In addition, a probability distribution 
function is associated to each photometric redshift. 

The primary template sets used are the four ob served 
spectra (Ell, Sbc, Scd, Irr) from lColeman et al.l fl980) com- 
plemented with two observed starburst sp ectral energy dis- 
tribution (SED) from lKinnev et all (1 19961 ). These templates 
have be en optimised using th e VVDS deep spectroscopic 
sample (|Le Fevre et al.l 120051 ). Next, an automatic zero- 
point calibration has been carried out using spectroscopic 
redshifts present Wl, W3 and W4 fields. For spectral type 
later than Sbc , a red dening E(B — V) = to 0.35 using the 
ICalzetti et all (2000) extinction law is applied. 

Although photometric redshifts were estimated to i < 
24, we selected obj ects brighter than i = 22.5 to limit out- 
liers (see Table 3 of lCoupon et al.|[2009h . Star-galaxy separa- 
tion is carried out using a joint selection taking into account 
the compactness of the object and the best-fitting templates. 

Finally, using the photometric redshift, the associated 
best-fitting template and the observed apparent magnitude 
in a given band, one can directly measure the fc-correction 
and the absolute magnitude in any rest-frame band. Since 
at high redshifts the fc-correction depends strongly on 
the galaxy SED, it is the main source of systematic er- 
rors in determining absolute magnitudes. To minimise the 
fc— correction uncertainties, we derive the rest-frame lumi- 
nosity at a wavelength A using the object's apparent magni- 
tude closest to A(l+z) according to the redshift of the galaxy 
(the procedure is described in Appendix A of lllbert et al.l 
2005). For this reason the bluest absolute magnitude selec- 
tion takes full advantage of the complete observed magni- 
tude set. However, as the u-band flux has potentially larger 
photometric errors, we use M 9 -band magnitudes. 



2.3 Sample selection 

We select all galaxies with i' < 22.5 outside masked re- 
gions. At this magnitude, our photometric redshift samples 
are 100% complete, for all galaxy spectral types, and our 
photometric redshift errors are small compared to our red- 
shift bin width in the redshift range < z < 1. 



CFHTLS Wide - all galaxies (I < 22,5) 



6 http : //terapix . iap . f r/article .php?id_article=841 

7 ftp : //f tpix . iap . f r/pub/CFHTLS-zphot-T0007/ 
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Figure 1. A two-dimensional histogram of the number density 
in the M g -z plane is plotted for all galaxies in the CFHTLS wide, 
T0007 release. The red line shows our adopted thresholds in ab- 
solute luminosity and redshift. 



Table 1. Number of galaxies per redshift bin outside the masked 
regions. 





Wl 


W2 


W3 


W4 


0.2 < z < 0.4 


36 307 


12 121 


23 426 


10 945 


0.4 < z < 0.6 


100 586 


31 586 


67 413 


32 778 


0.6 < z < 0.8 


165 215 


48 596 


105 020 


59 137 


0.8 < z < 1.0 


149 960 


45 783 


96 737 


52 423 



Here we chose M g < -20.64 (M B - 5 log h < -19.9) and 
divide this selection into four redshift bins: 0.2 < z < 0.4, 
0.4 < 2 < 0.6, 0.6 < z < 0.8 and 0.8 < z < 1.0. A large bin 
width (Az = 0.2) ensures a low bin-to-bin contamination 
from random errors, but still reduces substantially the mix- 
ing of physical scales at a given angular scale that occurs in 
an angular catalog without redshift information. 

The sample selection is illustrated in Figure [T] As a 
means of estimating the approximate completeness of our 
samples, we measured the mean absolute magnitude M g in 



each of the redshift bins and found M„ 



-21.4 ±0.1. With 



http : //www . cf ht .hawaii . edu/~ arnouts/lephare .html 



this selection we have more than 10,000 objects in each bin 
in each Wide patch. The detailed number of objects in each 
sample is summarized in Table [T] 



3 MEASURING GALAXY CLUSTERING 

3.1 Two-point angular correlation function 

As a consistency check with previous work, and as a means 
to verify our lower-order counts in cells moments, we mea- 
sure the two-point angular correlat i on fun ction w(Q) for our 
samples using the lLandv fc Szalavl l|l993l ) estimator, 
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w{9) 



n r (n r 



1) DD 



1 DR 



Tid(nd — 1) RR n d RR 



+ 1, 



(2) 



where, for a chosen bin from 9 to 9 + 59, DD is the number 
of galaxy pairs of the catalog in the bin, RR the number of 
pairs of a random sample in the same bin, and DR the num- 
ber of pairs in the bin between the catalog and the random 
sample, rid and n r are the number of galaxies and random 
objects respectively. 

A random catalog is generated for each sample with 
the same geometry as the data catalog using n T = 10 6 
which is at least 6 times and at the most 90 times na for 
a given bin and patch. We measure w in each field using 
the "Athe na" tree c odtPI in the angular range 9 £ [0.001, 1] 
degrees divided into 15 bins spaced logarithmically. 

In § 13. 41 we explain how the individual measurements of 
w(9) from the four different fields are combined. 



can use the following recur sion relation (|Szapudi fc Szalavl 
ll993bl ; [CoTombi et al'llioooh : 



71 — 1 



i\ (n - k) S n -k Fk 



n * — ' (n — k)\ k\ 

fe=i 



with 



£ = 6 = F 2 /F? - 1 



(8) 



and N c = iV£. 

In practice, one usually measures, for a cell of size 9, the 
count probability Pn{9), the probability for a cell to contain 
N galaxies, as described below. Once Pn(9) is measured, it 
is straightforward to calculate the factorial moments using 



F k = p nN 



(9) 



3.2 Higher-order moments 

Consider a distribution galaxy counts-in-cells which is a ran- 
dom sampling of a continuous underlying density field p(r) 
where r is the coordinate of the object. The factorial mo- 
ments of the number of galaxies in a set of random cells of 
size 6 are closely related to the moments of the field through 

F k = { iV lfel > = ( N{N - 1) • • • (N - k + 1) ) = N k p k (3) 

where { N ^ ) represents the k th factorial moment. Fi = 
(N) — N represents the average cell count and p — p ★ We 
is the underlying density field smoothed with the window 
function of the cells. 

The quantities of interest, however, are the hierarchical 
amplitudes Sk defined as 



Sk 



6 



e 2 



where 



1 



dn ... dr fc £ fc (n, 



(4) 



(5) 



is the fc-point correlation function averaged over the cell an- 
gular surface V = s or volume V — v according to whether 
statistics are considered in two or three dimensions, respec- 
tively. In particular, 



6 = S / d01 (102 10(01,02). 



(6) 



The parameters S3 and S4, proportional to the skew- 
ness and the kurtosis of the distribution of the numbers 
of galaxies in cells, quantify deviations from the Gaus- 
sian limit. On a more general level, the hierarchical am- 
plitudes Sk are expected to vary slowly with angular scale 
as a consequence of gravitational instability both in the 
weak ly non-linear regime, as predicted by perturbation the- 
ory l|juszkiewicz et all 1 19931 ; iBernardeau et aL 20021) and 
the hi ghly non-linear regime as shown by I Davis fc Peebles! 
lll977f) by assumin g stable clustering plus se lf-similarity (see 
also Ipeeblesl Il98d ; IBernardeau et all |2002| . and references 
therein) . 

To relate hierarchical moments to factorial moments one 



http: / /www2 . iap . f r/users/kilbinge/athena/ 



3.3 Measuring counts in cells 

The main issue in counts in cells measu rements is that large 
scales are dominated by edge effects (ISzapudi fc Colombil 
1996; IColombT et al. 1998) due to the fact that, as cells have 
finite size, galaxies near the survey edges or near a masked 
region have smaller statistical weights than galaxies away 
from the edges. To correct for the se defects, we use "BMW- 
P N " (Black Magic Weighted P w ; IColombi fc Szapudill2012l ; 
iBlaizot et ail 120061 ). 

This method is based on the fact that counts in cells do 
not depend significantly on cell shape for a locally Poissonian 
point distribution. This latter approximation allows one to 
distort the cells near the edges of the catalog and near the 
masks in order to reduce as much as possible edge effects 
discussed above. In practice, the data and the masks are 
pixellated on a very fine grid. In this case, one can consider 
all the possible square cells of all possible sizes. For each of 
these cells, an effective size is given which corresponds to 
the area that is contained in the cell after subtraction of the 
masked pixels. If the fraction of masked pixels in the cell 
exceeds some threshold /, the cell is discarded. Then, the 
center of mass for the unmasked part of the cell is calculated 
and the algorithm finds the corresponding pixel it falls into. 
As one cell per pixel is enough to extract all the statistical 
information, the code selects among all the cells of equal un- 
masked volume targeting the same pixel the most compact 
one, i.e. the one which has been the least masked. 

Note that the current implementation of BMW-Pjv 
works only in the small angle approximation regime. This 
is not an issue in our analyses as the largest field Wl has an 
angular extension smaller than 9 degrees. 

We used the following parameters in BMW- Pat: the 
sampling grid pixel size was set to 1.93 arcsec, which corre- 
sponds to a grid of size 16000 x 13944 for Wl. This very high 
resolution allowed us to probe in detail the tails of the count 
in cells probability at scales as small as m i n = 0.0011 deg. 
We measured counts in square cells of size 9 for 15 angular 
scales, logarithmically spaced in the range 9 £ [0 m in,0max] 
with max = 1 deg (see Tables [Bljl . 

The allowed fraction of masked pixels per cell / is a 
crucial parameter. If / is too large, the cells become too 
elongated, invalidating the local Poisson approximation. On 
other hand, if / is too small, edge effects become prominent. 
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We investigated the following set of values of / for Wl in 
the redshift bin 0.6 < 2 < 0.8: / — 0, which corresponds to 
a traditional counts in cells measurement with square cells 
not overlapping with the edges of the catalog or the masks, 
and / = 0.25,0.5,0.75. We found that / as large as 0.5 
is needed to ensure that edge effects are minimised in the 
scaling range probed by our measurements. Larger values of 
/ do not change the results significantly so we kept / = 0.5 
to reduce as much as possible cell shape anisotropies. 



3.4 Combining fields 

Factorial moments obtained from equation ([9]), where counts 
in cells are measured on a single field Wi, are unbiased esti- 
mators. However, concerning counts in cells, the quantities 
of interest are the averaged two-point correlation function 
£ and the hierarchical moments S n , which are non linear 
combinations of F n , as in equation s @ and © , and hence 
biased statistics (e.g. ISzapudi et al]|l999ft . To minimise the 
bias, it is preferable to first perform a weighted average of 
the factorial moments over the four fields with weights pro- 
portional to effective area S e s,i, 



F n = 



(10) 



where F ni i is the n th factorial moment measured on field i 
and S e s } i its effective area (equation [1]) - and then, to use 
equations (Q) and © to compute £ and S n from F n . 

We have decided to use the same weighting procedure 
as in equation (|10[) to combine the individual correlation 
functions Wi{9) measured in each field and derive a com- 
bined w(0). This simpler approach leads to slightly biased 
measurements but is sufficient for our purpose. 

To calculate error bars, we exploit the fact that we have 
at our disposal four independent fields to estimate statisti- 
cal uncertainties directly from the data, using the dispersion 
among the four fields. The calculation of such errors is po- 
tentially complicated, as each field has a different area. For 
the factorial moments, which are unbiased, the errors are 

(AF k ) 2 = -. y^Ses^Fkj-Fkf/y^See,*. (11) 

^fields - 1 ^ 



Details are provided in Appendix [A] 

Rigorously speaking, a full error propagation formula 
shoul d be used to compute the statistical u ncertainties of 
w(9 ) (ILandv fe Szafavlll993l ; iBernsteinl 1 1994T ) . and of £ and 
S n (|Szapudi et al.l [1999ft . However in our case an involved 
procedure such as this is unnecessary since the uncertain- 
ties on the estimated errors are very large due to the fact 
that we have access to only four independent fields. To keep 
the approach simple while preserving, in practice, sufficient 
accuracy, we compute the errors on w(9), £ and S„ as for 
the F k , by simply replacing F k ,i with Wi(Q), |j or Sk,i in 
equation pip (i is the field number) , where the quantities £i 
and Sk,i are derived from Fk,i using equations (0 and |J5J. 

The hierarchical moments and their associated error 
bars are presented for each redshift bin in Tables IB1I 



3.5 Transformation to physical co-ordinates 

Tr ansformation to physica l co- ordinates is treated follow- 
ing lGroth fc Peebles! ill977ft and lFrv fc Peebles! l | 1978ft . orig- 
inally due to iLimberl (| 1953T ) . Let n be the spatial density of 
galaxies and 4> De the selection function, that is, the prob- 
ability that a galaxy at comoving distance x or redshift z 
is included in the catalog. Summed over distances in a spa- 
tially flat universe, the projected number density M per solid 
angle in a given redshift range 2 m i n < z < 2 max is 



M— I a 3 x 2 dxn(p = 



n(z) dz, 



(12) 



where a photon at redshift z arrives from comoving distance 
dz' 



x(z) — cH 1 



o ^n m (i + z') 3 + q a 



(13) 



The density of galaxies per solid angle per redshift interval 
is 

r \ dN - 3 2 da; 
dzdil dz 

The factor na 3 is constant, and we take a(z = 0) = 1. Each 
sample has a characteristic depth weighted by geometry and 
redshift distribution, 



(14) 



(15) 



* _ J n{z)x{z)dz 
J n(z) dz ' 

where the integrals run from 2 m i n to z max . The angular den- 
sity scales as M ~ (a:*) 3 - 

The projected two-point function w(0) is related to the 
spatial two-point function as 



N 2 w{9) = / n(zi) dzi 71(22) dz2 £(2:12). 



(16) 



The correlation function is important only at small separa- 
tions, for which X12 is 



x 2 12 = (^) 2 Az 2 +x 2 A9 2 
y dz ' 



(17) 



Changing variables to z = ^(zi + Z2) and Az = 21 — 2:2, and 
defining y as y = (dx/dz)Az/x9 leads to 

rt 2 W (9) = J dzn 2 {z) J ^Lj(x9(l + y 2 )^ 2 ). (18) 

To the extent that x* is the characteristic scale for x, the 
second moment scales as w(9) = x*~ J Fct (9x*); that is, the 
comb ination x*w is a function of 9x* (cf. iGroth fc Peebles! 
1 1977ft . 

The projection depends on the form of the spatial corre- 
lation function. The observed correlation function is approx- 
imately a power law; including the stable clustering time 
dependence, £ = a 3 ~ J (x/xo)' J , we have 



w{9) = Cj 
where 



, — I 7 



J d2 n(z) 2 a 3 - 1 x 1 -' (dx^)" 1 



[jdzn(z)Y 

Ay _r(|)r(2ri)^ 



(l + l/ 2 )-T/ 2 



r(i) 



3.67909, 



(19) 



(20) 



approximated for 7 w 1.8. Thus, the angular two-point func- 
tion is also a power law, 
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Table 2. Projection parameters. Columns give the redshift 
range, scaling distance and n-point amplitude projection factors 
F p (n) (eq. [23]) for n = 3, 4, 5. 



z\ 


22 


x* (Mpc) 


*p(3) 


*i>(4) 


F P (5) 


0.2 


0.4 


1280 


1.09 


1.28 


1.61 


0.4 


0.6 


1910 


1.02 


1.06 


1.12 


0.6 


0.8 


2530 


1.02 


1.05 


1.10 


0.8 


1.0 


3010 


1.09 


1.25 


1.46 



){0) = TI xl 



~t a 1 — i 

11 



(21) 



with power index shifted by 1. 

The projected n-point function is given by the extension 
of equation (|16[) . 



ni dzi . . . n k dzfc 



(22) 



For hierarchical model where the n-point correlation func- 
tion can be expressed as a sum over products of 2-point 
correlation functions: £ n = ^ Qk (sum over connected 
terms), where Qk are the reduced amplitudes and £y = 
^2(xi, xj). After a change of variables to the average redshift 
and k — 1 differences, the angular correlation function is also 
found to be hierarchical, with amplitude qk = QkF p (k), 



F p (k) = 



"V 



h It 



n k dz 



3-7 1-7 \ fe-1 



/a" 1 'i 1 '\ 
V dx/dz / 



(23) 



(24) 



Table [2] shows numerical results for the redshift distribution 
of the data with 7 = 1.8. Variations from field to field are 
generally less than a percent. Changing 7 by 10% affects 
the results by at most a few percents, and considering the 
perturbation theory time dependence £(x) ~ D 2 (z) instead 
of a 3-7 by less than 1%. 



3.6 Effect of photometric errors 

tn order to quantify the effect of the photometric errors 
on the values of the S„, we investigate how important 
the variatio n is in the pro j ection factors. To do so we use 
results from lCoupon et al.l (|2009l ) who compute the photo- 
metric redshifts for the CFHTLS Wide T0004. They use a 
template-fitting method to compute photometric redshifts 
calibrated with a large catalog of 16,983 spectroscopic 
redshifts from the VVDS-F02, VVDS-F22, DEEP2 and 
the zCOSMOS surveys. Their method includes correction 
of systematic offsets, template adaptation, and the use of 
priors. Comparing with spectroscopic redshifts, they find a 
photometric dispersion in the wide fields of 0.037-0.039 and 
at i' < 22.5. We then take a dispersion Az = 0.038 (1 + z), 
and we get the numbers presented in Table 3 for x*, Fp2, 
Fps, Fpi, Fp$. We see that the errors on the photometric 
redshifts introduce a very small difference in the projected 
numbers and we can consider that they are negligible 
comparing to the error that we get from the field to field 
variance. 
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1.000 
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Figure 2. The two point angular correlation function, w(9) (sym- 
bols) compared with the measurement of Coupon et al. (dotted, 
dashed and dot-dashed black curves). The red curves represents 
the linear theory prediction taking into account linear bias, as 
discussed in Section l4.2.3l 



4 RESULTS 

4.1 Two-point angular correlation function 

As a first consistency check, Figure [5] shows our mea- 
surements of the two-point correlation function w(8) in 
four redshift b ins. W e compare three bins with those of 
ICoupon et al.l (|2012r i (dashed lines). On the scales shown 
in this Figure, agreement between the two measurements is 
excellent. This is to be expected: at i < 22.5 the precision 
of photometric redshifts in CFHTLS release T0007 is sim- 
ilar to the previous release, T0006, used in ICoupon et alj 
(2012). We over-plotted the dark matter predictions for the 
linear part of the two-point correlation function from which 
we determine the value of the bias in the different redshift 
bins. These are presented in Table 4. One can note that our 
values for the bias are close to unity do not strongly evolve 
with time. 

Considering the average two point correlation (defined 
by equation [6} on Figure [3] we see, as expected, the same 
behaviour as w{6) i.e that the clustering increases as the red- 
shift decreases. We can also use this measurement to verify 
that the Landy & Szalay estimator measurements are con- 
sistent with BMW-PN measurements. To do this, we fit w{9) 
by a simple power law w = A9 1 ' using a x 2 to estimate 
the amplitude A and the power 7. Our results are plotted 
on the left panel of Figure [3] 

Using the formula [5] with square cells, we derive an ap- 
proximate analytical expression for the average two point 
correlation function: 



4x2 (3-7)/2 
W ~ r Ad 

(3- 7 )(4- 7 ) 



(25) 



Comparing with a numerical integration, we found that 
this formula works for 7 in the range [1.6, 2.0] (our values of 
7 are 1.9, 1.8, 1.7 and 1.65 for the redshift bins 0.2 < z < 0.4, 
0.4 < z < 0.6 , 0.6 < z < 0.8 and 0.8 < z < 1.0 respectively) 
with a precision that increases with 7. The differences in w 
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21 z 2 x* (Mpc) F p (3) F p (4) F p (5) 



0.2 


0.4 


1272 


1.09 


1.29 


1.64 


0.4 


0.(5 


1905 


1.02 


1.07 


1.14 


0.6 


0.8 


2530 


1.02 


1.06 


1.11 


0.8 


1.0 


3005 


1.09 


1.25 


1.48 



Table 3. Projection parameters obtained using photometric redshifts and the dispersion in the estimated redshift of Az = 0.038 (1 + z). 
Columns give the redshift range, scaling distance and n-point amplitude projection factors F p (n) (eq. 1231 ) for n = 3, 4, 5. 



1000 
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0.001 0.010 0.100 

scale (degrees) 
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Figure 4. Hierarchical moments 53, 54 and S$ measured in the 
DR7 release of the SDSS for galaxies selected with < z < 0.4, 
18 < r < 21 and M r < —20.7 (dashed lines) compared to mea- 
surements in the CFHTLS (solid lines). 



between the two methods ranges from 8.5% for 7 = 1.6 to 
less than 1% for 7 = 2. The results are presented on the right 
panel of Figure [3] We can see that this prediction matches 
very well the measurements. 



4.2 Higher-order moments 

4.2.1 Comparison with SDSS DR7 

We now consider our measurements of the hierarchical mo- 
ments. As first consistency check, we compare our S n s mea- 
sured in the CFHTLS to measurements made from the sev- 
enth data release of the Sloan Digital Sky Survey (SDSS- 
DR7) using an i ndependently-written counts-in-cells code 
iRoss et all l|2007l ). This code was used with the default con- 
figuration parameters proposed by the authors. 

Since the SDSS and CFHTLS are by their nature very 
different surveys, we can only make this comparison for the 
lowest redshift bins in the CFHTLS. We therefore chose 
galaxies with < z < 0.4 18 < r < 21 and M r < -20.7 
(comparable to IRoss et "al]|2007t ) and applied the same se- 
lection to galaxies in the SDSS. Our results are plotted in 
Figure [4] where we show the hierarchical moments obtained 
from the four fields of the CFHTLS (solid lines) compared to 
measurements in the SDSS (dashed lines). For the CFHTLS, 
error bars are computed according to the method described 
in § 13.41 for the SDSS, error bars were computed using a 



jackknife method (e.g.. | Scranton et al.l I2OO2T ) . For the de- 
tailed method see lRoss et al. (|2006bl ). 

As a consequence of the extremely fine pixelisation grid 
we employ, we are able to measure our S n to much smaller 
angular scales than in the SDSS (such a measurement would 
be computationally demanding for a large-area survey like 
the SDSS). In the common angular range, the agreement be- 
tween the CFHTLS and the SDSS measurements is excellent 
except at the largest scales. In this latter regime, however, 
our measurements are expected to be less precise, due to 
finite size effects when approaching the angular size of the 
fields. 



4-2.2 Field-to-field variations 

Figure [5] shows the third-order cumulant S3 as a function 
of scale in different redshift bins for the four different fields. 
At low redshifts, z < 0.6, the dispersion between the four 
fields is low; however, above z ~ 0.7 W3 begins to be sys- 
tematically higher than the other fields at all scales. At the 
highest redshifts, the difference between W3 and the other 
fields is significant. The same behaviour is observed for the 
other cumulants. To illustrate furthermore this point Fig- 
ure compares the hierarchical moments as a function of 
scale for the various redshift bins for the cases when all the 
fields are combined together and when W3 is excluded. The 
presence of W3 is in fact crucial because when it is included 
there clearly seems to be a redshift evolution of the hierar- 
chical moments as indicated by the left panel of Figure [5l 
while removing W3 makes the redshift dependence insignif- 
icant as shown on the right panel. We now investigate this 
effect in more detail. 

The explanation of the behavior of the cumulants at 
the highest redshift bin can be understood from an exam- 
ination of Figure [7] which shows the counts in cells distri- 
bution function for the combined fields, W1+W2+W3+W4, 
for W1+W2+W4, and for W3 alone. Clearly, the presence of 
W3 dominates and perturbs considerably the high end tail of 
Pn{0). A careful inspection of the W3 field reveals that this 
is due to the presence of one or several rich clusters in W3 
comp ared to the three other fields (see e.g. IColombi et al.l 
119941 ). Due to the finiteness of the sample there is a value 
Amax(#), above which Pn(Q) becomes zero. This function 
Amax(#) obviously relates to the richest clusters of the sam- 
ple: for instance a cell of size 8 located in the densest part 
of one of the richest cluster of the sample should contain 
N ula , x (6) galaxies. Similarly, Pn(9) for TV close to A max is 
entirely determined for a given smoothing scale 9 by the 
density profile of a rich cluster of galaxies (cf. Figure 5 of 
IColombi et~ai1ll994l ). 

To illustrate furthermore this point, Figure [8] shows the 
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Figure 3. The symbols on the left panel represent the measurements of w(8) in the four redshift bins. The lines are the fits obtained 
using a \ 2 assuming a power law for the two point correlation function. The symbols on the right panel show the averaged two point 
correlation measured using counts in cells. The lines are the predicted values using formula (31) with the parameters of our w(6) fit. 



largest cluster in W3 at 0.8 < z < 1.0 and the effect on 
S3 of subtracting it from the field. As expected, removing 
the cluster greatly reduces the skewness over a large angular 
range. Indeed, moments of order k of the count probability 
are increasingly sensitive to the high-iV tail of Pn{0) with 
increasing k. Removing the largest cluster reduces the high 
N tail of Pn{9) and hence reduces the skewness. This hap- 
pens here only on a finite range of angular scales because, at 
smaller (larger) scales, one or several clusters with different 
concentration parameters from those of the one we selected 
dominate the high N tail of the count probability. In fact, 
detailed examination of the four fields in the redshift bin 
0.8 < 2 < 1 shows that there are four rich clusters in W3 
and one in Wl, which affect the cumulants and explains the 
discrepancy between W3 and the other fields. 



4-2.3 Significance of field-to-field variations 

To test if the discrepancy between W3 and the three other 
fields is statistically significant we decided to compute the 
theoretical expectation of the errors on £, 53 and Si in or- 
der to compare them to the errors estimated as described in 
§ 13.41 and as shown in Table IB1I This is made possible by 
the use of the package FO RCE (FORtran for Cosmic Errors, 
IColombi fc Szapudil 200 if) which e xploit s the analytical cal- 
culati ons of ISzapudi fc Colombil £1996) and ISzapudi et al.l 
(1999) to provide quantitative estimates of these statistical 
errors due to the finiteness of the sample, given a number of 
input parameters that specify a prior model. 

All the parameters needed for such a calculation were 
estimated from the survey self-consistently. These are the 
area S of the survey, the cell area s, the average galaxy 
number count N, the averaged two-point correlation func- 
tion over a cell £, the values of S n 's for 3 ^ n ^ 8 and the 
averaged two-point correlation function over the survey, 



this cluster has coordinates (216.5, 55.9) in radec J2000. 



Redshift bin 


b 


0.2 < z < 0.4 


1.17 


0.4 < z < 0.6 


1.21 


0.6 < z < 0.8 


1.27 


0.8 < 2 < 1.0 


1.37 



Table 4. Bias in the different redshift bins. 



£s= £j J d9 1 d9 2 w(9i,9 2 ). 



(26) 



Except for £s, all the parameters above can be extracted 
directly from the data, including S„ , which measured up to 
the eighth order. 

The calculation of £s is however more complex. To per- 
form it we use linear perturbation theory to derive a three- 
dimensional correlation function for a ACDM model. Then, 
we compute the angular two point correlation function in 
the various redshift bins using equation (1161) , from which we 
can derive numerically the integral (|26|) I 11 ! The result is mul- 
tiplied with the square of the bias between the galaxy and 
the dark matter distributions. To estimate this linear bias, 
we fit the measu red two-point correlation function using 
the "halo-model" (IScoccimarro et ai1l200ll: iMa fc Frvl[2000l : 
Peacock fc Smithl2000jJCoorav fc Shethl2002l ), exactly as in 
Coupon et all l|2012l )l 12 l The values are given in Table 14X31 

Figure [9] compares the theoretical errors on £, 5*3 and 
5*4 derived with FORCE to the measured ones using the dis- 
persion from the four fields. Three models are considered, as 
detailed in the caption of the figure, to reflect uncertainties 
on the theory. These theoretical predictions are expected 
to be valid only when relative errors are small compared to 



11 w{6) was computed using the halo-model module of C0SM0PMC 
(http://cosmopmc.info), while £5 was estimated by Monte- 
Carlo simulation using 10 6 pairs to perform the integral. 

12 using again C0SM0PMC. 
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Figure 5. The four panels represent measurements of 53 in the four fields of the CFHTLS from the lowest redshift bin 0.2 < z < 0.4 to 
the highest one 0.8 < z < 1.0. The red line represents the combined S3 computed using the method described in Section 13.41 We can see 
that at low redshift all the fields behave identically while at higher z, S3 in W3 is systematically higher over the whole scale range. 



unity as consequence of the propagation error technique used 
to derive it. Consequently, predictions at scales beyond 0.1 
degrees, where edge effects start to dominate in the theory 
are expected to be less accurate. With all these limitations 
in mind, one can see that the agreement between theory 
and measurements is excellent except in the redshift bin 
0.6 < z < 0.8 where the relative error seems to be slightly 
underestimated/overestimated by the measurements when 
comparing to the theory for £ and the hierarchical moments 
respectively. 

The remarkable agreement between theoretical errors 
derived using the standard ACDM model and the measured 
weighted dispersion between the four fields in the redshift 
bin 0.8 < z < 1 demonstrates that in fact W3 is not special, 
from the statistical point of view. This is confirmed by the 
examination of Figure 1101 which shows the hierarchical mo- 
ments as functions of redshift for two typical angular scales: 
8 = 0.012 deg, which corresponds in practice to the non- 
linear regime, and 6 — 0.22 deg, which is the largest possible 



reliable scale in order to probe the weakly non linear regime. 
On this plot, one sees that the S n measured after removing 
W3 are compatible within 2er with those obtained with all 
the fields. 



4-2.4 Evolution with redshift and comparison with 
perturbation theory 

We have demonstrated that W3 does not represent a "spe- 
cial event" and that measurements with and without W3 
are compatible within 2a. This demonstrates that the red- 
shift evolution trend in the non-linear regime observed on 
left panel of Figure |§] and on Figure [TU] is indicative but 
inconclusive; this statement holds when angular scales are 
converted to Mpc. From now on, we do not consider fur- 
thermore W3 separately and present results only for the full 
combination of the four fields. 

Figure [11] shows the hierarchical moments in three di- 
mensional space as functions of physical scale after applying 
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Figure 6. Combined Sn's as functions of angular scale in the four redshift bins when using W1+W2+W3+W4 together (left panel) and 
when W3 is excluded (right panel). Points which are not plotted correspond to negative values. 




Figure 7. Counts in cells distribution function Pjv(9) as a function of N in the last redshift bin 0.8 < z < 1.0. Three cases are considered: 
the combination W1+W2+W3+W4 using the weighting scheme as in equation (1 1 I t (thick grey curves), the combination W1+W2+W4 
(red dotted curves) and W3 alone (solid black curves). Each curve corresponds to a different scale; for clarity we keep only one scale out 
of two. 
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Figure 8. A rich cluster in W3. The upper left panel represents the distribution on the sky of galaxies in W3 for 0.8 < z < 1. The red 
square encloses the most prominent concentration from a visual inspection. The upper right panel illustrates the effect of masking this 
concentration on S3 for W3 in the highest redshift bin. The lower left panel shows a coloured image of the sky in the region delimited 
by the red square, which reveals a cluster of red galaxies. On the lower right panel, the galaxies of the same region are represented in 
the same field with symbol colour representing redshift. We use a wider redshift bin (0.7 < z < 1.0) to be sure not to cut galaxies that 
belong to the cluster due to photometric redshift errors. 



the transformation described in Section 13.51 The measure- 
ments show the same behavior as before i.e. a plateau at 
small scales corresponding to the highly non-linear regime, 
and then tend to another lower plateau at large scales cor- 
responding to the weakly non-linear regime. At scales ~ 10 
Mpc, we start to be sensitive to finite size effects as we ap- 
proach the size of the fields. 

We now compare the me asurements direct ly to 

perturbation theo ry predictions (lJuszkiewicz et al.l Il993l; 
Ir3ernardeaulll994al lbl; iBernardeau et a~ 20021 ) : 



34 



+ 7i, 



(27) 



oP t 60712 62 7 2 2 
54 = 132y + Y 7l + 3 7l + 3 72 ' 



„ PT 200575880 1847200 6940 2 

65 = „„^„,» — 1 — „„„„ 71 + -^rJi 



305613 3969 

235 3 1490 50 
+ — 71 + -^-72 + -77- 



(,:! 



10 



7i72 + ^73 



(28) 



(29) 



27 ' 63 

with 7 P = d p lncr 2 (i?)/d p ln_R, where a 2 (R) is the linear 
variance of the 3D matter density field smoothed with a 
spherical top-hat window of radius R. 

The agreement between perturbation theory predictions 
and the measurements at large scales is excellent and even 
better in a striking way, as it can be seen on in Figures [10] 
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Figure 9. Comparison of theoretical errors to measured ones. The different panels from the left to the right represent the relative error 
on respectively the averaged two-point correlation function. 53 and Sa comparing the different models available in FORCE (black curves) to 
the errors measured from the dispersio n over the four fields as ex plained in Se ction [ 3.41 (red curve). For the th eoretical predictions, three 
different models are available, SzS from | SzaPudi fc Szalavl dl993al lbt) . BeS from iBernardeau fc Schaefferl dl992l ^ and Per uses perturbation 
theory predictions of Bcrnardca 3 l|l994a| , ll99cfr Each row of panels represents a different redshift bin going from the lowest one on the 
top to the highest one on the bottom. Some curves are interrupted where the models fail and give negative error bars. 
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Figure 10. Evolution of S3, S4 and S5 with redshift at non-linear and a weakly non-linear angular scales. The black symbols linked 
with solid lines and those linked with dotted lines correspond to 9 = 0.012 deg for W1+W2+W3+W4 and without W3 respectively. 
The hashed grey region represent 2tr around the solid black curves. Similarly for 9 = 0.22 deg in red. In addition, perturbation theory 
predictions for the scale 9 = 0.22 deg are given in blue with linear bias (diamonds) and without (triangles). 





Figure 11. Three-dimensional Sn's in the four redshift bins, together on the same plot (upper left panel) and represented individually 
(three remaining panels). The dashed and dot-dashed black curves on each panel gives dark matter halo model prediction. The coloured 
lines represent pure perturbation theory predictions using equations (22]), J2EJ an d JUJ, in two cases: one without bias and one with the 
most extreme value of linear bias, b = 1.37, found in the last redshift bin. 



© 0000 RAS, MNRAS 000, 000-000 



14 M. Wolk et al. 



and 1111 when simply taking into account effects of linear 
bias: 

Sn ^ = fc2(n-2) ' ( 30 ) 

with the values of b measured in Section 14.2.31 and 
given in Table 14.2.31 With such an agreement well within 
the error bars, we did not fit higher order bias terms 
l|Frv fc Gaztanaga!li"993h ; these are clearly consistent with 
zero at the level of approximation of our measurements. 

In addition, on Figure 1111 the full halo-model predic- 
tions for t he da rk matter are given using the method of 
iFrv et all l|201ll ). What is interesting to notice here is that 
the halo model predictions and measurements are consistent 
concerning the location of the transition between the weakly 
non-linear regime and the highly non-linear regime. 



5 CONCLUSIONS 

We have measured the hierarchical moments of the galaxy 
distribution 5*3, 5*4 and S$ in a volume-limited sample of 
more than one million galaxies spanning the redshift range 
z ~ 1 to z ~ 2. Our measurements were made in the latest 
release of the CFHTLS wide and cover an effective area of 
133 deg 2 . 

We are able for the first time to make robust measure- 
ments of the evolution of the hierarchical moments as func- 
tions of redshift within a large range of angular scales. The 
survey consists of four independent fields, allowing us to es- 
timate the cosmic variance at each bin and angular scale. 

Our main results are as follows: 

(i) In the redshift bin 0.2 < z < 0.4 for a sample with the 
same cut in absolute magnitude, our results are in excellent 
agreement with the SDSS demonstrating that our survey is 
a fair sample of the Universe; 

(ii) This is confirmed by a comparison of the field to field 
scatter of the variance of counts in cells and the predicted 
hierarchi cal moments with the analy t ic "cosmic erro r s" pre - 
dicted bv lSzapudi fe Colombil <|l99fJ ): ISzapudi et all |l999); 

(iii) The hierarchical moments that we measure show two 
regimes separated by a transition region: first, a plateau 
at small scales (< l/i~ 1 Mpc) corresponding to the highly 
non-linear regime, and then the data tend to another lower 
plateau at large scales (> 7/i _1 Mpc) corresponding to the 
weakly non-linear regime; 

(iv) At non-linear scales, the amplitude of the hierarchical 
moments increase with redshift: however, a more accurate 
error analysis demonstrates that our results are also consis- 
tent with no redshift evolution within 2a; 

(v) In the weakly non-linear regime measurements are 
fully consistent with the perturbation theory predictions 
when a linear bias term (measured from our data) is 
taken into account. Hi gher-order bias terms corrections 
(|Frv fe Gaztanagalll993T j are negligible compared to linear 
bias given the precision of our measurements; 

(vi) The position of the transition between the non-linear 
and the weakly non-linear regime 1 < R < 10 ft -1 Mpc is 
fully compatible with halo model predictions. 

A robust conclusion concerning the redshift evolution of 
the hierarchical moments in the non-linear regime, and the 



corresponding implications for galaxy formation and evo- 
lution, is challenging. This is because the richest clusters in 
our sample can dramatically influence higher order measure- 
ments. A complete understanding of this effect will require 
a detailed modelling of the counts in cells probability using 
the halo model formalism and taking into account the halo 
mass of the most massive structure in our survey. This will 
be the subject of a subsequent work. 
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APPENDIX A: WEIGHTED STATISTICS 

In this Appendix, we describe the way we estimate errors of 
the combined quantities. 

Let Xi, i = 1, . . . , n, be independent variables, { XiXj ) = 
(xi) (xj) for i / j, with a common mean {xi) = /i and 
variances {{x t — /i) 2 } = of . Conceptually, the obser- 
vations of some quantity in independent fields. Let Wi be 
a set of weights we define and use to compute a weighted 
mean 



n 

— £" 



X = > WiXi, 

Wtot 

i=l 

and dispersion 

n 

y^wjjxj - x) 2 , 



(Al) 



(A2) 



with Wtot = E w i- We seek a relation between the dispersion 
among the samples, Ax, and the uncertainty in the mean, 
Ax. The expectation of the weighted mean is 



1 n 

{ x, ) = — Wi { Xi ) = n, 

Wtot 



and the expectation of the weighted variance is 
1 n 

{{Ax) 2 ) = > u>i {x 2 — 2xxi + x 2 ) . 

Wtot ^-^ 



(A3) 



(A4) 



The terms in the sum can be expressed in terms of fx and 
the of, 

/ 2 \ 2 . 2 

( Xi ) = n +(Ji, 



{ XXi ) = ^ Wj ( XiXj ) 

Wtot ^ 

i 

— \wi(p 2 + af) + Wj[i 

Wtot L ^ 



2 . WjCT; 

M H 

Wtot 



(A5) 



(A6) 



and 



WjWfc {XjXk , 



1 X ^ 2 2 



(A7) 



Adding these together, the terms in fj, cancel, leaving 

Wtot L 



2w % o 2 Ej w l a l 



Wtot 



1 \ ^ 2 1 \ ^ 2 2 

> Wiffi =— > Wi Oi . 

Wtot^ W?_.^ 



(A8) 



The variance in x is similarly computed to be 

{(Ax) 2 ) = {x 2 )-{x) 2 = -^J2 w > a >- (A9) 

W to t 

^ 

In several cases this leads to a simple relation between Ax 
and Ax. If all the Wj are equal then 



<(A*n = 



n-l 



and 

i 

For Wi = 1/a 2 

ii a \2 \ n-l 
{{Ax) ) 



(A10) 



(All) 



(A12) 
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and 

<(AS) 2 ) = ^-L-^. (A13) 

The weights Wi — 1/of are the minimum variance estima- 
tor, and choosing the effective area as weights is the best ap- 
proximation to this0 Therefore, we use equation (|A1|) with 
weights Wi = S c H,i as the average over fields, and adopt as 
our estimate of uncertainty 

((A ^ )= «^>. (A14) 



APPENDIX B: S N MEASUREMENTS 



1 at least in the Poisson approximation, when Xi = F\ measured 
in field i 
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0.2 < 2 < 0.4 0.4 < z < 0.6 
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0.6 < 2 < 0.8 
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0.8 < 2 < 1.0 
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Table Bl. Measured S n and their error bars in the bins 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8 and 0.8 < z < 1, respectively on 
the upper left, upper right, bottom left and bottom right table. On each table, the first column refers to the angular size of the square 
cells used to perform the measurements. The next columns display successively 53 and its error bars, S4 and its error bars and S5 and 
its error bars. The counts in cells method used to perform the measurements is described in § 13.31 The way the four fields are combined 
and the corresponding error bars are computed is detailed in 13.41 
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